How to calculate correlation functions of 
Heisenberg chains 

Rob Hagemans*, Jean-Sebastien Caux* and Jean Michel Maillet^ 

*Instituut voor Theoretische Fysica, Universiteit van Amsterdam, Amsterdam, the Netherlands 
^ Laboratoire de Physique, Ecole Normale Superieure de Lyon, Lyon, France 

Abstract. We describe a method for calculating dynamical spin-spin correlation functions in the 
isotropic and anisotropic antiferromagnetic Heisenberg models. Our method is able to produce 
results with high accuracy over the full parameter space. 
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INTRODUCTION 

One-dimensional quantum systems stand out because the effect of interactions is unusu- 
ally strong, which turns them into a very suitable arena for the study of strong corre- 
lations. Another striking and very useful feature of certain one-dimensional systems is 
the fact that many of their properties can be calculated exactly, using methods that are 
unavailable for higher-dimensional systems. For a long time, however, the reach of such 
calculations has been limited to static properties. Indeed, calculating dynamical proper- 
ties within the framework provided by exact methods has for a long time been regarded 
as hard, possibly even intractable. This is changing, however. With the advent of exact 
expressions for form factors in one-dimensional spin models, the possibility is opened 
up to calculate dynamical quantities. 

One such quantity that is of considerable interest is the dynamical structure factor. 
This quantity, defined for a general spin operator S" on the jth site of an A'^-site chain as 
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(where the subscript c indicates we are to take a connected correlator) owes part of its 
importance to its role in describing inelastic neutron-scattering experiments on quasi- 
one-dimensional crystals, where it is proportional to the intensity of scattered neutrons 
at the given momenta and energies. 

Neutron scattering is used to probe the magnetic structure of materials. In a quasi- 
one-dimensional crystal, the coupling between sites in one direction has a much higher 
energy scale than in the two others. If one does measurements at an energy scale that is 
intermediate between these two scales, the system behaves as a one-dimensional system. 
In many such materials, the local magnetic moments have two low-lying levels, meaning 
that the systems can effectively be described as spin-^ chains, bringing us into a field 



of research that has been of theorists' interest for quite a long time: that of Heisenberg 
chains. 

Heisenberg's spin model, nowadays also known as the XXX model, originally con- 
sisted of a lattice of spins-^ with isotropic nearest-neighbour couplings. In one dimen- 
sion, this model was solved by Hans Bethe, who presented a methodLi] to determine 
its eigenfunctions and their energies. A more general spin model was introduced by R. 
Orbach[2]. This model, which is also the one we shall be concerned with in this article, 
features an anisotropy in the spin-spin interaction in the z direction (and is hence known 
as the XXZ model), as well as an external field parallel to that anisotropy. It is given by 
the Hamiltonian 

// = 7 £ S)S)^, ^8)8]^, +^iS)8]^i -\)- hS) , (2) 

where 7 is the coupling strength, A is the anisotropy parameter, h the external magnetic 
field and 5j, 8^-, 8) are the usual spin-i operators acting on chain site j. 
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Bethe' s method — known as the Bethe Ansatz — has been employed to find static 
properties of many one-dimensional quantum systems. However, due to the intricate 
form in which the eigenstates are given by the Bethe Ansatz, calculating such quantities 
as norms, overlap functions, and correlation functions is not an easy task. 



CALCULATION SCHEME 
Bethe 's Equations 

The Hamiltonian (O conserves the spin component in the z direction (as well as the 
total spin). Therefore, the Hilbert space it acts on separates into disjunct subspaces with 
different expectation values of the magnetisation. Throughout this paper, we shall spec- 
ify the subspace by the number of spins pointing downward, M := — (tIj^iS^^. 

Throughout this paper, the reference state is defined to be the state with all spins point- 
ing upward, |0) := ®^=i IT) - The chain length N is always assumed to be even and the 
boundary is periodic. If the Hamiltonian is defined as above, the energy of the refer- 
ence state equals —JhN/2. From this reference state, all states with M up to can be 
reached with a Bethe Ansatz; to reach the other states, one should start from a reference 
state with all spins down. 

Bethe's solution gives the un-normalised wave functions | {A }) in terms of a set of M 
parameters {A} known as rapidities. They can be calculated from the Bethe equations, 
a set of M coupled equations 
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with 



(^)„ (A ) = sinh( A + niC /2) with ^ = arccos A for XXZ, (4) 

(j)„{X)=X+ ni/2 for XXX. (5) 

In general, equations for the XXX case A ^ 1 follow from those for the XXZ case by 
taking the limit ^ ^ with 7 := A/^ kept finite. 

To calculate the many solutions of these equations, one makes them explicit by taking 
the logarithm of these equations, thus introducing a set of M quantum numbers {/}, 
which are integer for odd M and half -integer for even M. 
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where 

= 2varctan[tanhA/(tannC/2)''] for XXZ; (7) 

0„ ( A ) = 2 arctan [2 A /n] for XXX. (8) 
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There is an exclusion principle at work: since wave functions vanish identically for 
coinciding rapidities, no two rapidities must be equal. Because of this, for real rapidities, 
all quantum numbers 7/ are distinct. In that case, these equations can easily be solved up 
to high numerical precision by iteration. 



The string hypothesis and its limitations 

However, roots of Bethe's equations need not be real numbers. Already in Bethe's 
original paper ^ it is conjectured that complex solutions will form structures (later 
dubbed strings) of roots with equal real parts and imaginary parts equally spaced and 
distributed symmetrically around the real axis or, in the XXZ case only, the axis in/ 2, 
as follows: 

Ka = K + hi^{nj + \-2a) + \in{\-Vj)+i5L for a G {l, , (10) 

where j G {1 . . .A^^} denotes the string type, a G {1 . . .M{\ enumerates all strings of 
type k. Uj is the length and vj G {1,-1} the parity of a string. The string centre A^ 

always lies on the real axis. In this definition we allowed for a deviation 5^^; the string 
hypothesis is the conjecture that complex solutions of the Bethe equations take this form 
with 5L = 0{e-''^) with c> 0. 

For A = 1 all string lengths are allowed and have positive parity; yet, in general, what 
lengths are allowed and whatparities are associated to them depends on the anisotropy 
in a rather involved way flSlaH]: if we call the number of string types A^,, then, for 
J G {1 . . . A(s}, we have lengths and parities given by 
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v; = (-1)LK-i)^H (11) 



where Ns = mi + I, i is chosen such that nii < j < mt^i, and the series {y), (m) are found 
from 



j-i = , Jo = 1 , yi = yi-2 + v/jz-i , 

mo = 0, mi=Y^Vk for/ G {!,...,/} , (12) 
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in which v,- > 1 are the partial quotients of i^/n given by 
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where we terminate the continued fraction after a suitable number / of partial quotients. 
Note that V/ > 2 to avoid an ambiguity in this definition. In the presence of strings of 
length two or higher, indeterminacies of the type 5/5 pop up in the Bethe equations. 
These, however, can be cancelled by taking a product of the Bethe equations for all 
roots constituting a string, leading to what is known as the Bethe-Takahashi equations; 
in logarithmic form, they are 
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with 
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These equations give the string centres in terms of a new set of mutually distinct 
quantum numbers 1^. Note that the Mf. must satisfy Y^kLi ^k^k = M. 
The energies of the states thus defined are 
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E = -jy : h(\N-M) for XXX. (17) 



and their momenta are 



2k 

^ = ;r£5,^,+iM, + — (mod2;r). (18) 



It is known [[ll,0,0,[i8] that the string hypothesis is not generally valid. However, 
after finding a numerical solution to the log-Bethe-Takahashi equations (fT4l) . one can 
explicitly check whether or not it is a solution of the original Bethe equations ^ to 



sufficient accuracy, by calculating the deviation of each string root. To first order, this 
deviation is given by 
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where 
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in which any factor of zero appearing in numerator or denominator is to be left out. 

There are a few classes of string solutions which we can discard without performing 
any calculation: no even-length string with rapidity zero can be a solution of the Bethe 
equations, as the left-hand side becomes non-finite; furthermore, if there is more than 
one odd-length string of the same parity with rapidity zero, the exclusion principle is 
violated. Rapidities at zero can easily be found: if the quantum numbers are spaced 
symmetrically around zero, a zero quantum number implies a zero rapidity. The wave 
functions in these cases must be obtained by a limiting procedure on the Bethe wave 
functions, which we do not carry out here. 



Determinant expressions 

Although Bethe eigenstates are hard to work with, exact expressions exist for the 
norm and overlap ifioll of states and for form factors iU, 12] (i.e. matrix elements 



of local spin operators). By inserting a complete set of eigenstates, the structure factor 
can be expressed as a sum over squares of form factors of the M-particle ground state 
|Gm) and these eigenstates, 

1 ^ r°° I 

S''"{q,(0):=- £ e'?0--/) / dt e"^' (s"j{t) S^JO) 

= 2n £ \{GM\S''g\a)f5{(0-(0a) (21) 

The ground state itself is excluded from the sum as we calculate the connected correlator 
{S''S")^:= {S"S") - {S") (S"). Note that for the form factors to be non-zero a must be 
in the right subspace: for the form factor, with equal number M of down spins as the 
ground state; for S^, the space with (M — 1) down spins. The form factor can then be 
calculated using a determinant expression. For the two-particle contribution at a fixed 
momentum, such methods were used in 

We shall calculate the full contribution coming from multi-particle intermediate states 
over the full Brillouin zone. In terms of the string rapidities, the determinant expressions 



are EM 
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and 
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with the matrices H, H , and P given by 

H-„ = H,h = K,b Pah = forbe{l...n-l} 
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with 
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where notations such as Hab are a shorthand for Hj^^ f^pi, and the values A^q, .^j are 
defined by extension of equation (fTOl) . It should be noted that although the definitions of 
elements of H and H are the same, when H is used the number of rapidities A is one 
less than the number of /i, and therefore we have not yet defined all of its columns. The 
final column ofH is given by 
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The norm [jsli] is given by ||{A}|| 
matrix defined as 



= =det4>({A}) with the reduced Gaudin 
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Enumeration of states 



With 2^ states in the Hilbert space, the form factor expansion by itself does not 
suffice to calculate the structure factor for anything but very small systems. What saves 
the day is that the number of states over which we have to sum can be vastly reduced 

First of all, as we already noted, we only need a single M subspace. Secondly, 
we can divide this subspace further into classes of states, which we shall call bases, 
distinguished by the number and character of excitations they contain. For instance, 
there is always a base containing only real rapidities, which we could denote (■). Then, 
for the XXZ chain, the first different type of particle is a rapidity with negative parity; 
if we have one of them we have base (■, 1). A base with three negative-parity rapidities 
and a two-string would then be (-,3, 1), and so on. 

For each such base, we can determine bounds for the quantum numbers {/}. We do 
this by letting the largest rapidity in the log-Bethe-Takahashi equation ([14)) go to infinity; 
since the order of rapidities equals the order of quantum numbers, all quantum numbers 
must be lower than or equal to the one that leads to rapidity infinity. Allowed quantum 
numbers / for string type j thus satisfy 
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where it is understood that all P are integer for odd Mj, half-integer for even Mj. Let us 
call the highest allowed quantum number li,; then, the number of available positions for 
string type j is 2/i + 1 . If the equality li, = 1^^^^^^ is satisfied, as happens in the XXX 
case, we have a rapidity that is formally infinite. These rapidities merit special attention 
which we shall bestow on them later on. 

First, let us turn our attention to the lowest string type: the real rapidities. Since the 
ground state forms a Fermi interval in terms of the quantum numbers, any occupied 
position outside this interval is a particle and any empty position is a hole. As for the 
other string types, every occupied position there counts as a particle. 

It is observed that the sum of form factors in a base decreases rapidly with the total 
number of particles of each type (particles, holes, and higher strings). Therefore, we 
order our states along these lines and sum over their form factors progressively to 



get an approximate value for the structure factor. We can check the accuracy of this 
approximation by comparing the total sum to the sum rules 
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In the isotropic case, at zero field, we can also check against the sum rule for the first 
frequency moment, 

£ |(GM|S^|a)|' = f(l-cos^)(i-ln2), (31) 

a^GM 

giving an extra check on each slice of fixed momentum q. 



Infinite rapidities 

In the isotropic model, from a given solution of the Bethe-Takahashi equations at M 
down spins , one can generate a cascade of solutions atM+l,M + 2,...A^/2by applying 
the descending operator Sq, without changing the existing rapidities. The action of the 
descending operator is equivalent to adding an extra rapidity at infinity: 

5o |{r})= |{r,±-}) (32) 

The momentum of a state is left unchanged by this modification and the energy is raised 
by hAM = h'm the presence of an external magnetic field h. The form factors of such a 
state are also easy to calculate: we have the relations 

|({Ai}|5, |{7,±-})|' = (33) 
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(M|5;|{r.±-}>r= '--l— - (34) 



as can be shown explicitly from the determinant expressions (|22|) . (1231) or by application 
of the Wigner-Eckart theorem, noting the rotational symmetry of the Hamiltonian and 
the vector nature of the spin operator [13]. 



RESULTS 

As an example of the quality of results we can achieve with this method, in Figure [U 
we show the transverse correlation function for an isotropic chain of N = 320 sites at 
M = N /A. More results for anisotropic chains can be found in fl^]. 




FIGURE 1. Transverse correlation function S-+{K, oj) at A = 1, A? = 320, and M = N/4 = 80. The 
energy (0 is shown in units of the coupUng In this calculation, the sum rule was satisfied to 98.6%. 
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